Nonadiabatic quantum pumping in mesoscopic nanostructures 
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We consider a nonadiabatic quantum pumping phenomena in a ballistic narrow constriction. The 
pumping is induced by a potential that has both spatial and temporal periodicity characterized by 
K and fi. In the zero frequency (Q. = 0) limit, the transmission through narrow constriction exhibits 
valley structures due to the opening up of energy gaps in the pumping region — a consequence of the 
K periodicity. These valley structures remain robust in the regime of finite fi, while their energies of 
occurrence are shifted by about hQ,/2. The direction of these energy shifts depend on the directions 
of both the phase- velocity of the pumping potential and the transmitting electrons. This frequency 
dependent feature of the valley structures gives rise to both the asymmetry in the transmission 
coefficients and the pumping current. An experimental setup is suggested for a possible observation 
of our nonadiabatic quantum pumping findings. 
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PACS numbers: 73.23.-b, 73. 23. Ad, 73.50.Mx, 72.10.-d 

The phenomenon of adiabatic quantum pumping in 
mesoscopic systems has attracted much attention since 
the first proposal by Thouless . The pumping phenom- 
ena refer to net transport of charges at zero bias. Mech- 
anisms giving rise to quantum pumping involves cychc 
deformations of more than one parameter in the system. 
The adiabatic deformation causes the transmission of a 
finite amount of charges in one deformation cycle. For 
the cases when this finite charge transfer is quantized 
the pumping could be of iniportance in establishing a 
standard of electric current [31 ■ The first experimental 
demonstration of the pumping of charges was reported 
by Switkes et al. . They applied sinusoidal voltages to 
two metal gates that define the shape of an open quantum 
dot. The phase difference between the two metal gates 
is an important adjustable parameter for the pumping of 
charges. Subsequent theoretical studies on this pumping 
mechanism in quantum dots have invoked a double oscil- 
lating barrier model 0, |^ . More recently, it was pointed 
out that quantum interference also plays an important 
role in the pumping of charges Thus far, most of the 
studies have concentrated on the adiabatic regime Wi. It 
is legitimate then to explore the nonadiabatic aspect of 
the pumping phenomena. 

In this work, our purposes are threefold: to treat the 
pumping phenomena nonadiabatically, to analyze in de- 
tail a novel pumping mechanism found in our results, and 
to propose an experimental setup for a possible realiza- 
tion of the mechanism. Towards these ends, we have im- 
plemented a generalized scattering-matrix method that 
allows us to solve the time-dependent Schrodinger equa- 
tion to effectively all orders in the pumping potential . 
Hence our results are not limited by the adiabatic ap- 
proximation. Furthermore, in analyzing our results more 
closely, we are able to identify a pumping mechanism 
that is associated with coherent inelastic scatterings of 
the transmitting electron in the pumping region. That 
the above quantum transitions play a decisive role in this 



pumping phenomena shows the nonadiabatic nature of 
the phenomena. We stress that this pumping phenom- 
ena do not require an asymmetry in the system config- 
uration. The phase velocity of the pumping potential 
alone is sufficient to cause the pumping. 

The pumping potential V{x, t) we invoked acts upon a 
certain region of a ballistic narrow channel. Apart from 
the overall spatial envelope of the potential, V{x, t) is 
periodic in time t, with a period T — 27r/r2, and pe- 
riodic along the longitudinal location x, with a period 
ip = 2tt/K. We assume a simple envelope profile for the 
potential. 



V{x, t) = Vo cos{Kx - nt)Q (L/2 - |x|) 



(1) 



expecting it to have captured the essential physics of 
the pumping mechanism. Within the pumping region, 
the pumping potential has a phase velocity v = Q/K, 
right-going in this case, and the longitudinal dimension 
it covered is L. As long as L ^ Lp, the detail form of 
the envelope profile — such as a less abrupt profile — 
should not change the pumping features found in this 
work. Meanwhile, similar potential form has been con- 
sidered by O. Entin-Wohlman et al. j9| in their study 
on the acoustoelectric effect in a finite-length ballistic 
quantum channel. The potential was generated from the 
surface acoustic wave, and was assumed to be significant 
only inside the channel while it was totally screened in 
the terminals. Even though the potential we consider in 
this work is of similar form to that of theirs, the physi- 
cal regimes of interest are different. They considered the 
regime that incoherent processes occur frequent enough 
in the channel to sustain a well defined local distribution 
of the electrons 9] . 

Our interest, however, is in the coherent regime: 
that electrons traversing the channel can maintain their 
phase coherence without encountering any incoherent 
processes. It turns out that the incoherent and the co- 
herent regimes exhibit different transport characteristics. 
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For instance, in the incoherent regime, the pumping fea- 
ture was found to be most significant only in the vicinity 
of the subband threshold 9], whereas in the coherent 
regime, we find significant pumping features in the con- 
ductance plateau regions. 

To simplify our presentation, wc choose the energy unit 
E* = Ep, the length unit a* — l/kp, the time unit 
t* — h/E*, and the frequency unit fi* = 1/t*. Here Ep 
and kp represent, respectively, the Fermi energy and the 
Fermi wave vector in the reservoirs. The Hamiltonian is 
of the form H = Hy+H^(t), where Hy = -9W+w^y^ 
contains a transverse confinement, leading to transverse 
subbands with energy levels e„ = (2n -I- l)ujy and a sub- 
band level spacing Ae = 2u!y. The time-dependent and 
the longitudinal part of the Hamiltonian Hx (i) is of the 
form Hx{t) = -d^/dx^ + V{x,t). 

At zero source-drain bias, the chemical potential /i is 
the same in the reservoirs that connect to the two ends 
of the constriction. The pumped current consists of con- 
tribution from all electrons, given by 

2e f°° 

/ = dEf{E - m) [T^{E) - T^m . (2) 

Here f{E — /i) is the Fermi distribution function and 
— e is the charge of an electron. In addition, and 
are, respectively, the total current transmission co- 
efficients for the right- and left-going incident electrons 
and for a given incident energy E. Contributions from 
all sidebands, denoted by m, and all subbands, denoted 
by n, have been taken into account, given by 

n m 

where the primed summation includes only propagating 
components of the transmitted electrons. The T^'l^^.^^ 
denote the transmission coefficients for the processes that 
an incident electron in subband n, energy i?, passes 
through the pumping region and ends up in energy 
E -\- mil on the other end of the constriction. These 
coefficients are solved by a generalized scattering-matrix 
method that we had established to treat time-dependent 
scattering potential nonperturbatively We had ap- 
plied this method to cases of longitudinally polarized 
field ^3 and transversely polarized field |l3|. Essen- 
tially, in this method, we segmented the scattering po- 
tential along the longitudinal direction into slices, solved 
the scattering matrix of each slice, and constructed the 
total scattering matrix |^ . The energy derivative of the 
pumped current dl/d/j, is given, in the low temperature 
limit, by 

BT 2p 

— = --[T_.(m)-T^(a*)]. (4) 

In our numerical examples, we present the dependen- 
cies of the current transmission coefficients and dl/dfi. 
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FIG. 1: The current transmission coefficient as a function of 
X for Q/As = 0.2, 0.6, and 1.0 in Figs. l(a)-l(c), respectively. 
Other parameters are K = 0.15, and L — 150. Solid curves 
and dashed curves are the transmissions of the right- and 
left-going incident electrons, respectively, and for Vb — 0.002. 
Dotted curves are for Vo — 0.0. 

The values of the parameters used are consistent with 
that for a typical high mobility two-dimensional elec- 
tron gas formed at a GaAs — Ala;Gai_2;As heterostruc- 
ture, with E* = 9 meV, a* = l/kp = 79.6 A, and 
n* = E*/h = 13.6 Trad/s. We also choose Uy = 
0.007 such that the subband energy-level spacing Ae — 
0.014 (~ 0.126 meV). For the pumping potential, we 
chose Vb — 0.002, K = 0.15, and the longitudinal range 
L — 150 such that KL /2tt ^ 3.58. The chemical poten- 
tial /i is replaced by 

A;+2' 

where the integral value of X stands for the number of 
propagating subbands in the constriction. The depen- 
dence of the current transmission coefficients and 
on X are depicted by the solid and the dashed curves in 
Fig. 1, respectively. The dotted curves are for Vb = 
and angular frequency equals 0.2Ae, 0.6Ae, and Ae, 
respectively, in Figs. l(a)-(c). 

The general features presented in Fig. 1 are summa- 
rized in the following. In the absence of the pumping 
potential, the transmission coefficients and ex- 
hibit the well understood quantized features, reflecting 
the number of propagating subbands in the constrictions. 



3 



That both of the transmission coefficients are given by 
the same dotted curve reflects the symmetry in the trans- 
mission of the electrons with respect to their incident di- 
rections. This symmetry, however, is broken by the intro- 
duction of pumping potentials to the constriction. Upon 
the action of the pumping potential, a valley structure is 
developed in each of the plateau regions, but the energies 
of occurrence for these valley structures are different for 

and T^. A valley structure in Fig. 1 is characterized 
by a substantial drop in the transmission coefficient, with 
maximum drop |ATi„ax| > 0.5, and a large energy width, 
with AXvaiiey ~ 0.4A£. Valley structures for occur 
at lower energies than that for on the same plateau. 
Interestingly, the separation in energy, or AX, between 
the minimums of the valley structures and belonging to 
the same plateau for and T^, increases with the an- 
gular frequency Q. of the pumping potential. In fact, our 
numerical results show that AX = fi/Ae. For the case 
when Vl = Ae, the separation becomes so large that the 

valley structures overlap with the valley struc- 
tures of the next higher subbands. These valley struc- 
tures and the asymmetry in the transmission coefficients 
are the key results in this work. 

We note in passing that there are small oscillatory 
structures in Fig. 1(a) in the energy range between the 
subband threshold and the valley structure. These har- 
monic structures are associated with multiple scatterings 
of the transmitting electrons in between the two edges of 
the pumping potential. 

Before wc present a physical explanation for the valley 
structures, we want to make one more remark: all the val- 
ley structures in Figs. 1(a) and 1(b) are almost identical 
in their profile regardless of their differences in the ener- 
gies of occurrence and the incident direction of the elec- 
trons. Furthermore, even though the valley structures 
in Fig. 1(c) arc truncated at the threshold energies, the 
remaining profile compares well with those in Figs. 1(a) 
and 1(b). These ^-independent features prompted us to 
check the = case, and we find indeed the same valley 
structure profiles except that the transmission symmetry 
is restored, with both and curves falling on top 
of one another. This f7 = result is not presented here, 
but we learn from this result that the energies of occur- 
rence for the finite Q. valley structures of and have 
been shifted by about Vl/2 and —51/2, respectively, from 
their = counterparts. The = valley structures 
are certainly associated with the formation of an energy 
band gap inside the region acted upon by the potential 
V{x). It is conceivable then that the valley structures in 
the finite Q. regime are caused by the establishment of a 
resonant coupling between some degenerate states via the 
pumping potential. To drive this point home and to ob- 
tain an analytic expression for the energies of occurrence 
of the valley structures, we propose a two-component ap- 
proximation for the wavefunction along the constriction. 



given by 



ex.p{ikx — iEt) 

-|-Cexp[i(fc - K)x - i{E - Q)t\. 



(5) 



The above wavefunction has taken up a near resonance 
approximation that describes the resonant coupling of 
an electron with photons. Substituting Eq. (5) into a ID 
Schrodinger equation 

[-d^/dx^ + Vocos{Kx-m)\ i^{x,t) = id/dtij{x,t), 

and retaining only terms of the form as in %l){x,t), we 
obtain 



E 



-h Efc + £k-K ± V (n -I- £fe 



-K 



(6) 

where E is the energy spectrum of an electron coupled 
with photons in the system, and ek = k"^ ■ 

The energy spectrum has an energy gap at the reso- 
nant coupling condition given by £k = £k-K ± fi- The 
plus, and minus, sign in the resonant coupling condi- 
tion is for positive, and negative, fc, respectively. This 
in turn corresponds to right- and left-going states. The 
ID kinetic energy Sk at the center of the energy gap 
is EGap = [(-^/2)(1 ± ^/K'^)Y, which can be approx- 
imated by K'^/A ±0/2 when O < /T^. In the case 
for a constriction, the energies of occurrence will be at 
-''^Gap = N + EGap/Ae. Thesc Xcap's correspond to the 
minimums of the valley structures in Figs. l(a)-(c). More 
specifically, in Fig. 1(a), Xoap = 1.51,2.51,3.51 for the 
right-going electrons, and Xq^p = 1-31, 2.31, 3.31 for the 
left- going electrons. The matches of these values with 
the minimums of the valley structures arc remarkable. 

Besides the energies of occurrence, Eq. (6) also 
provides us an estimate for the energy gap, given by 
AXcap = Vo/Ae = 0.143 for our case. This AXcap 
is less than AXvaiiey ~ 0.4, the width of the valley struc- 
tures in Fig. 1. We have tried other cases of longer L and 
find out that AXvaiioy decreases with increasing L. Its 
value becomes 0.2 aX L = 400. In addition, the maximum 
drop in the transmission coefficient |ATmax| for the valley 
structures is very close to unity when L = 400. Thus we 
conclude that the seeming discrepancy between AXcap 
and AXvaiiey in Fig. 1 is resulted from the finite L ef- 
fect. The detail profile of the transmission coefficients, 
however, is given by our nonperturbative approach. 

We have shown that the valley structures in the fi- 
nite ri regime arc resulted from the coherent inelastic 
scatterings in the pumping region. The O-dependence in 
the energies of occurrence Xcap reflects the breaking of 
the transmission symmetry by the phase velocity of the 
pumping potential. 

In Fig. 2, we present the effects of the transmission 
asymmetry on the quantum pumping by plotting dl / 
as a function of X. The parameters are the same as those 
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FIG. 2: The energy derivative of the pumped currents, in 
units of 2e/h, as a function of X for Sl/Ae = 0.2,0.6, and 
1.0 in Figs. 2(a)-2(c), respectively. Other parameters are the 
same as those in Fig. 1. Solid curve is for Vo = 0.002 and 
dotted curves are for Vo — 0.0. 



that has taken full advanta ge o f the recently developed 
finger gate technology [l^ Il4l|. The basic structure is 
a narrow constriction defined out of a two dimensional 
electron gas by a pair of split-gates (SG). On top of this 
SG, and separated vertically by an insulating layer of 
submicron thickness, are two interdigitated finger gates 
(FG) that are oriented transversely. Each of these FGs 
has a period Lp. The pumping potential could be gen- 
erated by subjecting the FGs to AC bias of the same 
frequency while maintaining a phase difference between 
the two FGs. The control of the phase difference in the 
bias voltage between two neighboring gates was success- 
fully demonstrated by Switkes et al [3j in the / ~ 10 
MHz region. More recently, phase-shifter operating in 
the / « 10 GHz region has been fabricated using a 
0.6-/im GaAs process |0. Hence we believe that the 
suggested experimental setup, though poses a stringent 
challenge to the experimentalists, is within reach of the 
present nanotechnology. 

In conclusion, we have proposed and have analyzed 
in detail a nonadiabatic quantum pumping mechanism. 
We have demonstrated the robustness of such pumping 
mechanism due to its resonant coupling nature. And 
we have proposed an experimental setup for the possible 
realization of such mechanism. 

C.S.T. thanks Professor Y. C. Lee for his enthusiastic 
encouragement. This work was supported by NSC of the 
ROC under Contract No. NSC89-2112-M-236-001 and 
NSC90-2112-M-009-003. 



in Fig. 1. Major features shown in Figs. 2(a), 2(b) are al- 
ternate valley and broad peak structures. These are due 
to suppression in the transmission for electrons incident 
from the right, and left, electrode, respectively. The val- 
ley and the broad peak structures barely resemble each 
other in Fig. 2(a). However, the corresponding struc- 
tures in Fig. 2(b) are almost inversion images of one 
another. Since the valley structures for the two trans- 
mission coefficients overlap in Fig. 1(a), their contribu- 
tions to dl/dfji suffer from partial cancellation between 
themselves. This cancellation hampers the resemblance 
between the valley and the broad peak sturctures in Fig. 
2(a). On the other hand, the valley structures are well 
separated in Fig. 1(b), thus preserving the resemblance 
between the valley and the broad peak structures in Fig. 
2(b). The case of exact cancellations is presented in Fig. 
2(c), when the valley structures of overlap with the 
valley structures of of the next subbands. Only the 
first valley structure survives and the net flow of the elec- 
trons is in the same direction as the phase velocity of the 
pumping potential. Finally, we stress that the pumping 
effect is most significant in the conductance plateau re- 
gions. This is in sharp contrast to the pumping effect in 
the incoherent regime 0. 

For the possible realization of the pumping mechanism 
proposed in this work, we suggest an experimental setup 
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